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Abstract: Using a criterion proposed by Salpeter and standard solid-state 
physics, we have determined the Debye temperature of a solid in equilibrium 
with the electron gas surrounding it. The results obtained can have astro- 
physical applications in the determination of parameters of interstellar and 
interplanetary clouds. 
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Introduction 



Phase transitions of various kinds have been, and still are, actively con- 
tributing to important changes occurring in the world around us. The aim 
of the present paper is a contribution to the study of the conditions for the 
occurrence of the so-called "plasma-solid (PS)" transition. 

Simple physical reasoning shows that the PS transition must certainly 
occur in various astrophysical settings, such as in proto-planetary and proto- 
stellar clouds. On the one hand, it is widely known that the universe contains 
different kinds of plasmas. On the other hand, it is equally widely known that 
solid objects also exist in the universe. This implies that there must exist a 
physical region where a transition between the two regimes takes place. Note 
also that investigations of the PS transition could be interesting in non- 
astrophysical plasma physics, such as in considerations of a laser produced 
plasma in front of a metal target (cf. Dimitrijevic and Konjevic,1980). 

In a preliminary study of the PS transition (Celebonovic and Dappen, 
2000), we have determined the conditions for its occurrence using two simple 
idealized systems: a pure Fermi-Dirac and a pure Bose-Einstein gas. The 
object of the present calculation is again a determination of the conditions for 
the occurrence of the PS transition, but this time in a more realistic model. 
Our starting point is the criterion for the occurrence of the PS transition 
(which was not given that name) proposed by Salpeter (Salpeter,1961). 

Salpeter's paper was devoted to a thorough discussion of a zero-tempera- 
ture plasma. In the third part of that paper, he considered a system of 
positive ions of given charge and mass, rigidly fixed in the nodes of a perfect 
crystal lattice, and went on to estimate the zero point energy of the ions. He 
showed that the behavior of such a system can be described by the ratio 



where E Z:P denotes the zero point energy of the ions and Eq is the Coulomb 
energy. According to the analysis of (Salpeter, 1961) a PS transition occurs 
for / = 1. Calculations reported in that paper were estimates, and strictly 
valid only for T = K. In the calculations reported in the following,we have 
somewhat reformulated this criterion; namely, we have compared the energies 
per particle of a solid and the electron gas. 
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This means that in our formulation,the behaviour of the system is gov- 
erned by the ratio 

f m = ^ (2) 

where E PjS denotes the energy per particle of a solid, and E PjB is the energy 
per particle of the electron gas. Imposing the condition f m — 1 defines a 
line (or a region) in the phase space of the system where the energies per 
particle in the two phases are equal. In standard phase transition theory one 
of the conditions needed for phase equilibrium in a multi-phase system is the 
equality of the chemical potentials of its phases (for example, Landau and 
Lifchitz,1976). Salpeter's criterion is in fact equivalent to that statement, but 
under one supplementary condition - that the energy is independent on the 
particle number. Moving off this line (or out of the region) defined by f m — 1 
implies the occurence of a phase transition. Our calculations go beyond the 
original assumptions of Salpeter in at least two aspects: (i) we have not used 
estimates but exact calculations, and (ii) our results take into account the 
influence of temperature. 

Calculations 

Finding a general expression for the energy of a real solid that takes into 
account most (if not all) of its characteristics is a formidable problem in solid 
state physics (for example Born and Huang ,1968 or Davydov 1980). The 
complexity of the problem is caused by several factors:the frequency distri- 
bution of the oscillations of particles around their equilibrium positions, their 
mutual interactions and their large number (which is of the order of 10 23 ). 
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Various approximations to the complete solution exist; one of them is the 
Debye model of a solid. Within the Debye model, the energy per mole of a 
solid is given by the following expression (Born and Huang, 1968) 



E = N { u (v) + 9nk B T 
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where N denotes the number of elementary cells in a mole of the material, 
n the number of particles per elementary cell, u(v ) = jf = the static lattice 
energy per cell, 6 the Debye temperature, ks the Boltzmann constant, and 
in the following we will use the convention ks = 1. The second integral in 
Eq. (3) can be solved as (Abramowitz and Stegun,1972) 
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under the condition n >z 1, \x\ -< 2ir. These conditions are in (Abramowitz 
and Stegun, 1972) introduced for purely mathematical reasons. Physically, 
\x\ -< 2n implies T > ^.This will be used later in this paper,to express 
the Debye temperature of a solid in equilibrium with a plasma,as a function 
of the number density. Here, the symbol B 2 k denotes Bernoulli's numbers. 
Taking the static lattice energy per cell as the zero point of the energy scale, 
it follows from Eq. (3) that the energy per particle of a solid within the Debye 
model is given by 



E p , s = 9T 



where 
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Using Eq. (4) in Eq.(6), one finally gets 
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which implies that 
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Expanding explicitly it follows that 
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The energy per particle of a Fermi-Dirac gas of particles of number density 
n mass m and spin s is given by (Landau and Lifchitz,1976) 

3/2 

E >< = ¥^ F ^ M (10) 

All the symbols in this equation have their usual meaning. In particular, g = 
2s + 1, and F 3 / 2 (/?//) is the k = 3/2 case of a Fermi integral 

Fu{^)= / ^ -7^ vf ■ (11) 



1 + exp [{3 (e - //)] 

Analytical developements of this type of integrals into power series are known 
in the literature (such as Celebonovic,1998). The ratio f m as defined in eq.(2) 
can be formed from eqs.(8) and (10). After some algebra one gets the following 
general expression for this ratio: 

00 75 / n \ 2& 

fm = ( tfv»n/^ + 3 g (a ^ ( » ) ] (i2) 

This result has the interest of being mathematically general (within a partic- 
ular model of a solid). However, due to this generality hardly any physically 
useful conclusions can be drawn from it. Limiting the sum in eq.(12) to 
terms containing (^) 2 and inserting a known power series developement for 
the Fermi integral leads to the following approximation for the ratio f m : 

9n 20T 2 + 9 2 
Jm 2KT (12A 2 n 4 /3 + 57r 2 T 2) ^ 6 > 

The constant A in the preceeding equation occurs because of the expression 
for the chemical potential of the electron gas at T = K : /i = An 2 / 3 , and it 
is given by A = (37r 2 ) 2 ^ 3 h 2 /2m.The symbol K denotes the following combi- 

3 /2 

nation of constants: K = 2 i™^ 2h A ■ 
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Imposing the condition f m — 1 on Eq.(13), it becomes possible to deter- 
mine the conditions for the occurence of a PS transition. Solving Eq. (13) 
for the Debye temperature 6 gives 



= \^[2KT/ny[5n 2 T 2 + 12A 2 n 4/3 y[l 



90nT 



] (14) 



(5tt 2 T 2 + 12A 2 n 4 /3) 



We have thus obtained the equation of state of a system undergoing a PS 
transition. This equation links the relevant parameters of a plasma (number 
density and temperature) with those of the solid (Debye temperature) which 
is in phase equilibrium with it. 



Apart from being interesting from the point of view of pure statistical physics, 
this result can find astrophysical applications in studies of dust and gas 
clouds. Recent examples of observational studies of such clouds are for 
example (Jessop and Ward-Thomson, 2000) and (Grun,Kriiger and Land- 
graf, 2000). The particle number density and the temperature of an astrophys- 
ical plasma can be determined from observation using laboratory methods, 
such as the analysis of spectral line broadening (e.g.Dimitrijevic,1996) , or 
in some cases, in-situ measurements from space probes. Inserting these values 
into eq.(14) gives the possibility to determine the Debye temperature of a 
solid which is in equilibrium with the plasma surrounding it. A real physical 
example of such a system could be a condensing proto-planetary cloud, such 
as the one around the star /3 Pictoris. 

Consider Eq.(14) in the limiting case T — > O.Astrophysicallythis limit 
corresponds to an isolated low temperature interstellar cloud,or to a proto- 
planetary disk sufficiently distant from the central star. Mathematicallythis 
means that all functions of T in eq.(14) can be developed in series with T as 
a small parameter. Keeping only the first terms in these developements,one 
finally gets the following expression for the Debye temperature: 



Discussion 
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All the constants occuring in eq.(15) have been previously defined. 
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Inserting the values of all the constants which occur in Eq. (15), and re- 
taining only the biggest term,one finally arrives at 

9 = 67387.5T 5/2 n- 3/2 (16) 

This expression is the final approximation derived in this paper for the De- 
bye's temperature of a solid in equilibrium with the plasma surrounding 
it. Its applicability is limited by the condition ^ -< 2rr, which is built-in in 
eq. (4). Introducing this condition in eq.(16) leads to: 

T -< 0.00205621 x n (17) 

This is the physical limit of applicability of eq.(16).It is physically expectable, 
because of the limitations inherent to Debye's model on which the calculation 
leading to eq.(16) is based. Expression (16) has two possible applications. 
It can be used for the calculation of 9 if the temperature and the number 
density of the cloud are known. Recent observation shows that the tem- 
perature in the cold protostellar clouds in our galaxy is 10 < T[K] < 20 
(Cesarskyet al.,2000). Typical values of the density are in the interval 10 4 < 
n(cm~ 3 ) < 10 5 . Plotting the behaviour of eq.(16) within these limits shows 
that 9 < 1600 K. In a similar way, the pair of values n = 1.25 x 10 4 cm~ 3 , 
T = 15 K would lead to 9 = 630 K. Such a value is only slightly higher than 
the experimental value for the element Si. Note that this last result could be 
useful for the interpretation of observations, because emission from silicate 
grains has indeed been detected (Cesarsky, et al.,2000) Turning the argu- 
ment around, if solid particles are observed in a cloud, and if their chemical 
composition can be determined, it becomes possible to calculate the Debye 
temperature from the principles of solid state physics. If, in addition, tem- 
perature is known from spectroscopy, one can determine the value of the 
chemical potential, and finally, the number density of a cloud from Eq. (16). 
Further work along these lines is in progress and will be discussed elsewhere. 
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